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ABSTRACT 


The  purpose  of  this  thesis  is  to  lay  the  foundation  for  analyzing  structural 
coupling  effects  for  a  proposed  trolley  interface  between  a  ship  and  a  roll-on  roll-off 
discharge  facility  (RRDF).  Such  a  facility  could  allow  heavy  cargo  transfer  at  higher  sea 
states.  Previous  studies  have  analyzed  motions  assuming  that  there  is  no  structural 
coupling  between  the  trolley  and  the  RRDF.  A  mathematical  model  that  incorporates 
structural  coupling  is  developed  using  the  principle  of  virtual  work.  In  order  to  assess  the 
degree  of  necessity  for  the  proposed  model  we  conduct  a  systematic  series  of  numerical 
experiments.  In  these  calculations  we  model  the  trolley  through  a  generalized  stiffness 
coefficient  and  assess  its  influence  on  RRDF  motions.  It  is  shown  that  modeling  of 
structural  coupling  may  be  necessary  depending  on  the  relative  order  of  magnitude  of 
trolley  structural  rigidity  and  trolley  placement. 
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I.  INTRODUCTION 


One  of  the  salient  features  of  the  United  States  Military  approach  to  21st  Century 
Joint  Warfare  is  Maneuver  Warfare.  Maneuver  warfare  seeks  to  attack  adversary 
vulnerabilities  from  a  position  of  advantage  through  the  synchronized  application  of 
movement  and  fires.  This  is  one  of  the  reasons  for  the  creation  of  the  Naval 
Expeditionary  Strike  Group  (NESG).  Current  NESG  comprises  of  several  classes  of 
ships  such  as  the  LMSR,  CAPE-D,  CAPE-H  and  CAPE-T  to  deliver  the  necessary 
military  equipment  and  stores  to  support  its  mission.  In  situation  where  a  friendly  port  is 
available  close  to  the  objective,  these  ships  can  unload  their  payload  through  the  port 
facilities.  However  there  are  situations  where  a  friendly  port  is  not  in  close  proximity 
resulting  in  a  need  to  perfonn  at  sea  transfer  of  the  payload.  This  involves  ship-to-ship 
transfer  of  payload  from  these  large  ships  to  smaller  crafts  that  subsequently  delivers  the 
payload  to  shore. 

Current  ship-to-ship  transfer  involves  a  ramp  interface  for  heavy  equipment  such 
as  the  M1A1  main  battle  tanks.  This  technique  is  limited  to  operations  at  sea-state  3  due 
to  excessive  torsional  stresses  on  the  ramp  when  two  main  battle  tanks  are  located  in  the 
proximity  of  the  middle  of  the  ramp.  A  previous  study  on  a  new  method  of  transferring 
cargo  based  on  the  use  of  a  trolley-rail  interface  shows  potential  advantages  in  the 
reduction  of  torsional  stresses  and  weight.  This  thesis  lays  the  foundation  for  analyzing 
structural  coupling  effects  for  the  proposed  trolley  interface  between  a  ship  and  a  roll-on 
roll-off  discharge  facility  (RRDF),  and  determines  the  necessity  for  modeling  of 
structural  coupling  depending  on  the  relative  order  of  magnitude  of  trolley  structural 
rigidity  and  trolley  placement. 
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II.  BACKGROUND  INFORMATION 


Amphibious  operations  conducted  by  the  US  military  involves  large  amount  of 
equipment  and  stores  that  needs  to  be  unloaded  near  to  the  objective.  These  equipment 
and  stores  are  usually  kept  in  prepositioned  ships  at  various  parts  of  the  world  such  that 
the  response  time  of  the  forces  can  be  greatly  reduced.  However,  the  current  mode  of 
operations  requires  these  prepositioned  ships  to  unload  their  cargo  either  at  nearby  port 
facilities,  which  is  the  preferred  option,  or  in  the  event  that  such  ports  are  not  available, 
go  right  up  to  within  2  nautical  mile  of  the  beach  to  unload  their  cargo  onto  barges  or 
small  crafts  that  ferry  these  equipment  and  stores  ashore.  There  are  two  main  problems 
with  the  current  mode  of  operation.  First  of  all,  the  number  of  available  ports  for  the 
prepositioned  ships  to  unload  their  cargo  is  severely  limited.  Secondly,  the  number  of 
beaches  that  would  enable  these  large  prepositioned  ships  to  go  right  near  and  unload 
their  cargo  on  the  beaches  is  also  very  limited.  The  reason  for  this  is  that  the  transfer  of 
cargo  from  the  preposition  ship  to  the  barge  can  only  be  conducted  at  sea-state  three  and 
below  to  prevent  excessive  torsional  stress  on  the  ships’  ramp  when  unloading  heavy 
equipment  such  as  the  68-tonnes  M1A1  main  battle  tank  [Ref  1].  These  problems, 
coupled  with  the  fact  that  such  operations  results  in  extended  time  when  the  amphibious 
force  needs  to  defend  the  logistic  operations  meant  that  a  better  way  has  to  be  developed. 
This  is  the  main  driver  for  the  new  Trolley-Rail  system  proposed  by  Advanced  Design 
Consulting  Inc,  which  aims  to  enable  the  transferring  of  cargo  between  ships  at  higher 
sea  states.  Studies  on  the  response  of  the  proposed  system  under  seaway  has  yield 
positive  results  indicating  that  the  system  should  be  able  to  achieve  the  required  transfer 
rate  at  higher  sea  states  [Ref  2], 
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III.  MODELING 


A.  TROLLEY-RAIL  MODEL 

The  Trolley-Rail  system  proposed  by  Advance  Design  Consulting  Inc.  essentially 
comprise  of  two  parallel  rails  that  runs  from  one  ship  to  the  other  ship.  The  trolley  (or 
platform)  traverses  along  the  two  rails  via  attached  wheels  to  transport  equipment  and 
stores  from  the  prepositioned  ships  to  shuttle  ships  that  ferry  these  cargo  ashore.  For  the 
model,  the  other  ship  is  the  Roll-on  Roll-off  Discharge  Facility  (RRDF).  The  RRDF  is 
made  up  of  a  number  of  interconnected  modules  which  serves  as  the  platform  to  support 
the  rails  running  from  the  ship,  and  to  facilitate  the  movement  of  the  equipment  and 
stores  to  the  shuttle  ships.  A  schematic  of  the  Trolley-Rail  system  is  given  shown  in 
Figure  (1)  while  the  three-dimensional  model  is  shown  in  Figure  (2),  and  the  RRDF  is 
shown  in  Figure  (3): 


Figure  1.  Schematic  of  Trolley-Rail  System 
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Figure  2.  3-D  model  of  Trolley-Rail  System 


Figure  3.  Schematic  of  Roll-on  Roll-off  Discharge  Facility  (RRDF) 
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The  rails  are  each  modeled  as  a  beam  structure.  This  is  a  good  model  as  the 
length  of  each  rail  is  very  much  greater  than  its  width  and  thickness.  The  rails  are 
assumed  to  be  pinned  at  the  ship’s  end  while  the  other  end  is  free.  Each  rail  is  connected 
by  a  spring,  dashpot  and  actuator  to  the  RRDF  to  simulate  potential  active  control 
mechanism.  The  RRDF  is  connected  by  a  spring  and  dashpot  to  the  ground  to  simulate 
the  motion  of  the  seaway  [Ref  3].  The  trolley  will  effect  structural  coupling  on  the  two 
rails  and  this  is  being  modeled  as  a  spring  and  a  dashpot.  Two  lump  masses  are  also 
attached  at  the  trolley  interface  with  the  two  rails  to  account  for  the  cargo  loading  on  the 
trolley.  A  schematic  of  the  model  is  shown  in  Figure  3. 


Figure  4.  Model  of  Trolley-Rail  System 
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The  parameters  for  the  model  are  given  as  follows: 

L  =  length  of  each  Rail 
/  =  length  between  supports  along  Rail 
lp  =  length  of  trolley 

p  =  distance  along  the  beam  from  ship  end  to  the  center  line  of  the  trolley 

Md  -  mass  contribution  of  RRDF  at  the  support  location  of  the  free  end  of  beam 

k/i  =  stiffness  ratio  between  Rail  #  1  and  the  RRDF 

d  i  =  proportional  damping  ratio  between  Rail  #1  and  the  RRDF 

kn  =  stiffness  ratio  between  the  RRDF  and  the  ground  at  Rail  #  1  locality 

C12  =  proportional  damping  ratio  between  the  RRDF  and  the  ground  at  Rail  #  1 
locality 

kjj  =  stiffness  ratio  between  Rail  #  2  and  the  RRDF 

C21  =  proportional  damping  ratio  between  Rail  #2  and  the  RRDF 

k22  =  stiffness  ratio  between  the  RRDF  and  the  ground  at  Rail  #  2  locality 

C22  =  proportional  damping  ratio  between  the  RRDF  and  the  ground  at  Rail  #  2 
locality 

vi(x,t)  =  transverse  deflection  of  Rail  #1 

Z](t)  =  vertical  displacement  of  RRDF  at  locality  of  Rail  #1 

zsi(t)  =  vertical  displacement  of  the  sea  wave  at  locality  of  Rail  #1 

V2(x,t)  =  transverse  deflection  of  Rail  #2 

Z2(t)  =  vertical  displacement  of  RRDF  at  locality  of  Rail  #2 

zS2(t)  =  vertical  displacement  of  the  sea  wave  at  locality  of  Rail  #2 
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B.  DERIVATION  OF  MATHEMATICAL  MODEL 
1.  Equation  of  Motion 

Having  developed  the  schematic  of  the  model,  the  next  step  is  to  derive  the 
equations  of  motion  for  the  model.  Lagrange  equation  was  employed  for  this  purpose. 
First  of  all,  the  kinetic  energy  of  the  system  is  given  by 
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Working  on  the  premise  of  small  angles,  the  angular  velocity  of  the  RRDF  may  be 
approximated  by, 
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where  w  is  the  width  of  the  RRDF. 


Hence,  the  kinetic  energy  is  given  by 
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The  strain  energy  of  the  system  is  given  by 
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For  the  rails  which  are  modeled  as  beam  structures,  the  assumed  mode  method  is 
employed  to  approximate  the  transverse  vibration  of  each  beam,  given  by 

N 

v,  (x,t)  =  YJVl(x)aM) 

i= 1 
N 

v2(x,t)  =  YJVi(x)pi(t) 

;=i 

where  t//l  (x)  is  the  assumed  modes  that  satisfies  the  geometric  boundary  conditions 
of  the  beam  and  N  is  the  number  of  modes  that  was  used  in  the  computation. 

To  reduce  the  mathematical  complexity,  indicial  notations  are  utilized  to  represent  the 
transverse  deflection  of  the  beam.  Therefore  the  transverse  deflection  of  each  beam  is 
given  as 

v^x,t)  =  y/iai 
v2(x,t)  =  i//ipi 

where  the  subscript  i  represents  a  summation  from  1  to  N. 

Hence,  the  kinetic  energy  of  the  system  is  represented  in  the  form, 

T  =  pd{v/.(x)«,(0}  {y/j^djit)}  dx 

+  j £  pA  {<//, (x)fli (/)} { (//. (x)Pj (/)} dx 

+\Mp  \{¥Sp)^M)){¥j(p)^St))\  +  ^ MP  [{i^Cp)A(0}  {¥j{p)Pj  (0} 
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while  the  potential  energy  of  the  system  is  represented  in  the  form, 


U  =  ^El{y/'Xx)ai(t)}{y/'i(x)aj(t))dx  +  ^ku  (t)- z,f  +^knzx2 

+ 1 J,'  El { y/ (x)fi (/)} { y/j  (x)/l (/)}  dx  +  ^k2l  [(//, (I)fi (t)  -  z2  f  + | £22z22 

+  [vt (P)a, (0  -  Wi( P)P, (Of 


Applying  the  Lagrange  equation, 


d_ 

dt 


f  dT^ 


dq, 


■M-a=o 

dq,  dqi 


the  generalized  coordinates  for  the  model  are  chosen  as, 


q  \z\  zi  ^2  •  •  ■  aN  A  A  •  •  •  PN 


The  various  differentials  in  the  Lagrange  equation  are  given  by 

'  dT 


d_ 

dt 


V5ii  J 


i(z,-z;) 

/  W  J 


d_ 

dt 


dT 


-Md+\jd  1(2, -2,) 

v/  w  ) 


d_ 

dt 

d_ 

dt 


^  pT  i\T  7T  J  V-T  S—  j  ' 


1 

2" 
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8T  _8T  _8T  _8T  _  Q 
<9zj  dz2  8aj  8j3; 

^-  =  ~kn[  Vi  OK  (0  -  ]  +  knzi 

8z  j 

=  (kn+kn)zl-kny/i(l)al(t ) 

8U 

—  =  (*2l  +  ^22  )  Z2  -  k2l  ViWPiit) 
dz2 

|^={r 

+  ^M  [i//  (x)a,(0- ](//,(/) 

+  kP  [Vi (P)a, (0 - ^ GO. Pj (0] ^(/O 

^ = a:  ei  w>  w]  <**} 

+  ^2i[|//,(x)/?,(0-^2]'//(0 
-  kp  [  V,  (p)a,  (0  -  Wj  (P)Pj  (0]  Vj  ip) 

The  principle  of  virtual  work  is  used  to  determine  the  generalized  force  The 
RRDF  and  the  beams  are  assumed  to  have  very  small  movements  as  follows: 

<5zp  c5z2,  Svx(/,t),  and  Sv2(/,t)  at  location  jc  =  / 

Svx(p,t),  and  Sv2(p,t)  at  location x  = 
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The  virtual  work  done  due  to  dampers  resisting  force  is  then  given  by 
8Wd  =  -c12zxSzx  —  cxx  [ij  -v1(/,t)][(5’z1  -8vx(l,t)\ 

-  c22z2Sz2  -  c21  [z2  -  v2  (/, 0]  [Sz2  -  8v2  (/,  0] 

-  CP  [T  ( P ,  0  -  v2  ( p ,  0]  [Svx  ( p ,  t)-Sv2  ( p ,  0] 


SWd  =[-(c12  +cnK  +cui//l(l)di(t)\Sz] 

+  [-(c22+c2X)z2  +c2xVfi(l)Mt)]*z2 
+  [c,  xzxy/j  (/)  -  cnyrt  {l)y/j  (l)dj  (t)]  Sat 


+ 


Vi  (p)«,  (0  -  Vj  ( P)Pj  (0}  Vj  (p) 


Sa, 


+  \c2^iVi  (0  -  c2i^-  (0^-  (0  A  co] 


+ 


{ V ,  (pH  (0  -  Vj  ( P)pj  (0}  Vj  (p)] 


Hence,  the  generalize  forces  due  to  the  dampers  are  given  by 
Qi  =-{cn+cn)zl+cnvl(l)ai(t) 

Qz2  —  —  (C22  ~*~C2l)^2  C2\Vi  (0A(0 

Qt,  =  ci  (0^i  +  [-Ci  (Ovo  (0 + (p)vo  (p)]  « /  (0  -  (p)vo  (p)]  A  (0 

0A  =  C21  Vi (l)z2  - \cpV,(p)Vj (P)] aj (t)  +  [~c2]y/l (/)^7 (/)  +  cp^(. (p)y/j (p)\ pj ( t ) 
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For  the  external  forces,  the  virtual  work  done  is  given  by 
SWf  =  -M . 


(z  —  z  ^ 

*sl  ^ s2 

' 8 7  -8z2 

-Jjd 

(  8zx  -  Sz2  N 

l  2  J 

{  2  J 

l  w  J 

Jo  pA 7^i(0^(x,o^ 


10 '  l 

X 

jo'  7 


*l  x 

■Jo  pA  —  zs2(t)Sv2(x,t)dx 


SWJ  =  - 


n 


i 


7  Md  +—Jd 

\4  w  j 


^  ^  d  +  _2  7/ 


W  ) 

r  ^  x 

J0  pd-zjt)if/i(x)dx 


(4-^2)^z  i 


i*  Z, 

j0  PAjZs2(t)Yi(x)dx 


Sa: 


Sfit 


Hence,  the  generalized  forces  due  to  external  forces  are  given  by 


\\pAX-Vi{x)dx 


770 


Qa=- 


\LQ  pAj<//i(x)dx 


z,2(0 
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Finally,  the  equation  of  motion  in  generalized  coordinates  fonn  is  given  as 


M{«}+[c]{«} +[*]{«} ={e} 


whereby  the  matrices  are  given  by 


M= 


Mr  M  r  0(2x2;y) 

-Mr  Mr 


J(2Nx2) 


M 


a(NxN) 


M 


, S(NxN) ) 


where  MR  = 


f\  1  ^ 

~lMd  +~Jd 

v  4  w  J 


/V  /V  pL  \ 

M*J  =J0  (X)dX  +  p¥iiP)¥j(P) 


(^11  ~*~^12  ) 

0 

(O(lxJV) 

0  ^ 

V(1  xN) 

0 

n  Ki ) 

O(lxjV) 

-ki\Vi(J)(\xN) 

(O(JVxl) 

^(iVxl) 

K 

a(NxN) 

Kb 

ap(NxN) 

0(Vx  1) 

-^21^/(0(jVxl) 

Kb 

ap(NxN) 

K„ 

p(NxN )  J 

/v  p  L 

where  Kaj  =  jQ  El y/"(x)y/'(x)dx  +  k, X\pi{l)\pj (/)  +  kpi//i(p)y/j(p) 
P a.  =  J,'  EIij/'(x)if/'j (x)dx  +  klxy/i  (l)y//  (/)  +  kpy/i  (p)y/j  (p) 

KPi]  =-kpVi(p)Vj(p) 
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[C] 


f 


V 


\ci2  +  cn) 

0 

cn¥i(l\uN) 

O(lxV) 

0 

—  (c22  +c2l) 

O(lxV) 

C2\¥id\lxN) 

cn¥i(0(Nxi) 

0(2Vx  1) 

ca 

“(JVxJV) 

^aP(NxN) 

O(TVxl) 

C2\¥lO(Ny.Y) 

^aP(NxN) 

S' 

X 

r  ^ 

where  Ca  =  -c, ,(//,  {l)y/j  (/)  +  Cpy/t  (p)i//J  (p) 
Caij  =  ~c2ly/,  (l)¥j  (/)  +  cpy/t  (p)y/J  (p) 

Cap,  =-Cp¥i(P)¥j(P) 


-MR&I-Zsl) 
M  R  ( z.s i  —  ^si ) 


C^'vl 


(JVxl) 


Qi 


.S' 2 


'(JVxl) 


1 


where  Qtj  =  -£  pA  —  y/i{x)dx 


2.  Assumed  Modes  for  Beam 

The  assumed  modes,  which  satisfy  the  geometric  boundary  conditions,  need  to  be 
used  to  calculate  the  generalized  mass,  damping,  stiffness  and  force  matrices.  The 
general  solution  of  the  beam  is  given  by  [Ref  4] 

y  =  A  cosh  fix  +  B  sinh  fix  +  C  cos  fix  +  D  sin  fix 
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The  boundary  conditions  are: 


At  x  =  0  < 


At  x  =  L  < 


rd v 

v  dx  jx=0 


V  dx~  J  x=L 


v  =  0 

rf2f=o 

.  dx 

(M  =  0) 

K'=» 

dx 

(M  =  0) 

o 

II 

(V  =  0) 

=^+c=o 

=>  ^  =  -c 

=  /?(^-c)  =  o 

=>  ^  =  c 

A  =  C  =  0 


=  p2  [B  sinh  ( pV)  -  D  sin  ( PL ))  =  0 


sin  (PL) 

B  = - P—PD 

sinh(/?T) 

=  p2  (B  cosh  (PL)  -  D  cos  {pL))  =  0 


fd2y) 

V  dx'  J x=L 


cos  (PL) 

=> B  = - yP-JzD 

cosh  (/?T) 

sin(/?T)  cos(/7L) 


sinh  ( PL )  cosh  ( PL ) 
=>  tan(/?Z)  =  tanh(/?T) 
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In  addition,  the  Eigenfunctions  y/i  utilized  in  the  assumed  mode  method  are 
orthogonal  to  each  other,  and  therefore  the  following  simplifications  can  be  effected  on 
the  matrices 

0  (i*j) 

pAL+{-M  p\//?(p)  (; i  =  j ) 

~  f  — ,  V/d/)'//,  (/)  +  kpi//i(p)i//j(p)  (i*j) 

“iJ  [EIAfL  +  kx  ly/i  (/)  +  y/i  (p)y/J  (p)  (i  =  j ) 

g  -k^y/^y/jil)  +  kpy/l{p)\//j(p)  (i*j) 

a,i  \ EI KL  +  k2 ,  <//,  ( l)i//j  ( /)  +  i//i  (p)y/J  (p)  ( i  =  j ) 
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IV.  NUMERICAL  EXPERIMENTS 


A.  METHODOLOGY 

Now  that  we  have  developed  a  mathematical  model  suitable  for  structural 
coupling  of  the  trolley/RRDF  interface,  we  need  to  develop  some  guidelines  with  regards 
to  its  applicability.  One  main  question  that  we  need  to  answer  is  whether  such  a  model  is 
indeed  necessary  for  this  problem.  In  order  to  assess  this  we  proceed  as  follows: 

1.  First  we  evaluate  the  motions  of  the  trolley/RRDF  assuming  that  there  is  no 
structural  coupling.  This  is  achieved  in  the  frequency  domain  utilizing  results 
obtained  by  WAMIT.  WAMIT  is  a  software  package  that  solves  the 
hydrodynamic  radiation/diffraction  problem  by  a  panel  discretization  of  the 
body.  It  is  a  3-dimensional  program  and  as  such  it  can  model  hydrodynamic 
interactions  between  multiple  bodies.  In  our  case  the  two  bodies  are  the  ship 
and  the  RRDF.  Since  we  allow  for  hydrodynamic  interactions  (in  other  words 
the  motion  of  each  body  depends  on  velocities  and  accelerations  of  the 
adjacent  body)  the  equations  of  motion  possess  12  degrees  of  freedom.  Six 
degrees  of  freedom  (surge,  sway,  yaw,  heave,  roll,  pitch)  come  from  each 
body. 

2.  Based  on  the  computed  body  motions,  we  evaluate  the  relative  motions  at  the 
interface  between  the  two  bodies;  i.e.,  at  the  trolley  end-points.  These  relative 
motions  are  subtracted  (taking  into  consideration  their  relative  phase  angle)  in 
order  to  obtain  the  angular  displacement  of  the  trolley  as  well  as  its  twist. 

3.  The  next  step  is  to  model  the  structural  rigidity  of  the  trolley  through  a 
generalized  stiffness  tenn.  The  relative  motions  of  the  trolley  are  then  used  to 
generate  the  reaction  force  in  the  trolley  due  to  its  stiffness.  This  force  is  then 
added  to  the  hydrodynamic  excitation  force  for  the  RRDF  and  the  process 
continues.  It  should  be  mentioned  that  such  coupling  procedure  is  not 
necessary  for  the  ship  since,  due  to  its  high  displacement,  it  is  unlikely  that  it 
will  be  appreciably  affected. 
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B.  RESULTS 

Results  are  presented  in  the  next  chapter.  All  results  are  presented  in  standard 
polar  plots  in  terms  of  the  RMS  value  of  the  indicated  response  (trolley  angle  or  twist). 
The  radial  coordinate  of  the  graphs  is  the  significant  wave  height,  while  the  angular 
coordinate  is  the  wave  heading.  Fully  developed  Pierson-Moskowitz  and  developing  or 
decaying  Bretchneider  wave  spectra  are  used  for  the  calculations.  Although  the  results 
have  been  obtained  for  long-crested  seas,  it  is  expected  that  the  main  conclusions  will  be 
valid  for  short-crested  seas  as  well.  The  trolley  stiffness  that  was  introduced  in  the 
previous  section  is  applied  towards  the  roll  motion  of  the  RRDF.  This  is  typical  for  a 
stern  configuration.  Similar  results  will  apply  for  other  trolley  placements  as  well.  For 
the  numerical  results  that  follow,  the  roll  stiffness  for  the  trolley  has  been  nonnalized 
with  respect  to  the  roll  stiffness  of  the  RRDF.  Naturally,  the  latter  is  a  function  of  the 
RRDF  displacement  and  its  metacentric  height. 
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V.  RESULTS 


A.  TROLLEY  ANGLE  AND  TWIST  SUBJECTED  TO  PIERSON- 

MOSKOWITZ  WAVE  SPECTRUM 

Figure  (5)  to  (10)  shows  the  RMS  value  of  the  trolley  angle  and  trolley  twist  when 
subjected  to  the  fully  developed  Pierson-Moskowitz  wave  spectrum  for  the  cases  where 
the  normalized  trolley  stiffness  was  set  to  a  value  of  0,  1  and  2.  180  degrees  corresponds 
to  head  seas  while  0  degrees  corresponds  to  following  seas.  The  anti-symmetry  in  the 
fore  and  aft  responds  is  due  to  the  placement  of  the  trolley  on  the  aft  of  the  ship.  The 
slight  anti-symmetry  in  the  Port  and  Starboard  responds  is  a  result  of  the  slight  anti¬ 
symmetry  of  the  shape  of  the  RRDF.  It  is  observed  that  the  trolley  angle  is  not  being 
affected  by  the  existence  of  the  trolley  stiffness  or  even  if  the  trolley  stiffness  is  doubled 
in  strength.  On  the  other  hand,  it  is  observed  that  the  twist  experienced  by  the  trolley  is 
significantly  affected  by  the  introduction  of  the  trolley  stiffness  coupling  tenn.  The  effect 
on  the  trolley  twist  is  amplified  further  when  the  normalized  trolley  stiffness  is  doubled. 


RMS  value  of  Trolley  Angle  (deg) 
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270 


0 


PM  -  Roll  Stiffness  =  0 


Figure  5.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Pierce-Moskowitz  spectrum  and 

zero  Roll  stiffness 
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RMS  value  of  Trolley  Angle  (deg) 
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Figure  6.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Pierce-Moskowitz  spectrum  and 

a  Normalized  Roll  stiffness  of  1 
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RMS  value  of  Trolley  Angle  (deg) 
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Figure  7.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Pierce-Moskowitz  spectrum  and 

a  Normalized  Roll  stiffness  of  2 
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RMS  value  of  Trolley  Twist  (deg) 


Figure  8.  RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Pierce-Moskowitz 

spectrum  and  zero  Roll  stiffness 
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RMS  value  of  Trolley  Twist  (deg) 
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Figure  9.  RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Pierce-Moskowitz 

spectrum  and  a  Nonnalized  Roll  stiffness  of  1 
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RMS  value  of  Trolley  Twist  (deg) 
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PM  -  Roll  Stiffness  =  2 


Figure  10.  RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Pierce-Moskowitz 

spectrum  and  a  Nonnalized  Roll  stiffness  of  2 
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B.  TROLLEY  ANGLE  AND  TWIST  SUBJECTED  TO  BRETSCHNEIDER 
WAVE  SPECTRUM 

Figure  (11)  to  (22)  shows  the  RMS  value  of  the  trolley  angle  and  trolley  twist 
when  subjected  to  the  Bretchneider  wave  spectrum  for  wave  periods  of  6  seconds  and  12 
seconds.  The  results  are  also  presented  for  the  three  cases  where  the  normalized  trolley 
stiffness  is  set  to  a  value  of  0,  1  and  2.  In  this  instance,  we  observe  a  similar  trend  in  the 
results  with  the  case  under  the  Pierson-Moskowitz  wave  spectrum  -  the  trolley  angle 
does  not  seems  to  be  affected  by  the  existence  of  the  trolley  stiffness  coupling  term  while 
the  trolley  twist  is  substantially  affected  by  the  introduction  of  the  trolley  stiffness  and  the 
effects  are  amplified  with  an  increased  in  the  strength  of  the  coupling  term. 
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Figure  11.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Bretchneider  spectrum  with  a 

wave  period  of  6  seconds  and  zero  Roll  stiffness 
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Figure  12.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Bretchneider  spectrum  with  a 
wave  period  of  6  seconds  and  a  Normalized  Roll  stiffness  of  1 
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RMS  value  of  Trolley  Angle  (deg) 
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Figure  13.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Bretchneider  spectrum  with  a 
wave  period  of  6  seconds  and  a  Normalized  Roll  stiffness  of  2 
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RMS  value  of  Trolley  Twist  (deg) 
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Figure  14.  RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Bretchneider  spectrum 
with  a  wave  period  of  6  seconds  and  zero  Roll  stiffness 
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Figure  15. 


B6  -  Roll  Stiffness  =  1 

RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Bretchneider  spectrum 
with  a  wave  period  of  6  seconds  and  a  Normalized  Roll  stiffness  of  1 
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Figure  16. 


RMS  value  of  Trolley  Twist  (deg) 
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B6  -  Roll  Stiffness  =  2 
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RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Bretchneider  spectrum 
with  a  wave  period  of  6  seconds  and  a  Normalized  Roll  stiffness  of  2 
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RMS  value  of  Trolley  Angle  (deg) 
90  30 
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B12  -  Roll  Stiffness  =  0 

Figure  17.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Bretchneider  spectrum  with  a 

wave  period  of  12  seconds  and  zero  Roll  stiffness 
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RMS  value  of  Trolley  Angle  (deg) 


B12  -  Roll  Stiffness  =  1 

Figure  18.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Bretchneider  spectrum  with  a 
wave  period  of  12  seconds  and  a  Normalized  Roll  stiffness  of  1 
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RMS  value  of  Trolley  Angle  (deg) 
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B12  -  Roll  Stiffness  =  2 

Figure  19.  RMS  Value  of  Trolley  Angle  (degrees)  subject  to  Bretchneider  spectrum  with  a 
wave  period  of  12  seconds  and  a  Normalized  Roll  stiffness  of  2 
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B12  -  Roll  Stiffness  =  0 


Figure  20.  RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Bretchneider  spectrum 
with  a  wave  period  of  12  seconds  and  zero  Roll  stiffness 
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RMS  value  of  Trolley  Twist  (deg) 
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Figure  2 1 .  RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Bretchneider  spectrum 
with  a  wave  period  of  12  seconds  and  a  Normalized  Roll  stiffness  of  1 
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RMS  value  of  Trolley  Twist  (deg) 
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Figure  22.  RMS  Value  of  Trolley  Twist  Angle  (degrees)  subject  to  Bretchneider  spectrum 
with  a  wave  period  of  12  seconds  and  a  Normalized  Roll  stiffness  of  2 
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c. 


INFLUENCE  OF  TROLLEY  STIFFNESS  ON  TROLLEY  TWIST 


The  observations  from  the  previous  sections  are  better  illustrated  through  Figure 
(23)  to  (28),  which  shows  the  comparison  of  the  responds  of  the  trolley  angle  and  trolley 
twist,  with  and  without  the  trolley  stiffness  coupling  term.  The  almost  symmetrical 
response  of  Port  and  Starboard  direction  is  expected  with  the  slight  variation  being 
accounted  for  by  the  slight  anti-symmetry  of  the  shape  of  the  RRDF.  It  can  clearly  be 
seen  that  there  is  substantial  differences  in  the  responds  of  the  trolley  twist  when  the 
stiffness  coupling  tenn  is  introduced.  This  responds  is  magnified  when  the  coupling  term 
is  doubled.  In  addition,  it  is  observed  that  the  quartering  seas  results  in  higher  trolley 
twist. 
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Figure  23.  Difference  in  average  vertical  trolley  twist  angle  for  a  Normalized  Roll  stiffness 
of  0  and  1  (subjected  to  Pierce-Moskowitz  wave  spectrum) 


90 


30 


39 


PM  -  Roll  Stiffness  -  0  vs  2 

Figure  24.  Difference  in  average  vertical  trolley  twist  angle  for  a  Normalized  Roll  stiffness 
of  0  and  2  (subjected  to  Pierce-Moskowitz  wave  spectrum) 
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Diff  RMS  Trolley  Twist  (deg) 
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B6  -  Roll  Stiffness  -  0  vs  1 

Figure  25.  Difference  in  average  vertical  trolley  twist  angle  for  a  Normalized  Roll  stiffness 
of  0  and  1  (subjected  to  Bretchneider  spectrum  with  a  wave  period  of  6  seconds) 
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Figure  26.  Difference  in  average  vertical  trolley  twist  angle  for  a  Normalized  Roll  stiffness 
of  0  and  2  (subjected  to  Bretchneider  spectrum  with  a  wave  period  of  6  seconds) 
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Diff  RMS  Trolley  Twist  (deg) 
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Figure  27.  Difference  in  average  vertical  trolley  twist  angle  for  a  Normalized  Roll  stiffness 
of  0  and  1  (subjected  to  Bretchneider  spectrum  with  a  wave  period  of  12  seconds) 
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Diff  RMS  Trolley  Twist  (deg) 
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Figure  28.  Difference  in  average  vertical  trolley  twist  angle  for  a  Normalized  Roll  stiffness 
of  0  and  2  (subjected  to  Bretchneider  spectrum  with  a  wave  period  of  12  seconds) 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  thesis  has  laid  the  foundation  for  the  analyzing  of  structural  coupling  effects 
for  the  proposed  trolley  interface  between  a  ship  and  a  Roll-on  Roll-off  Discharge 
Facility  (RRDF)  through  the  development  of  a  mathematical  model  utilizing  Assumed 
Modes  and  Virtual  Work  techniques,  which  incorporates  the  structural  coupling  effects  of 
the  trolley  to  the  RRDF.  Through  a  systematic  series  of  numerical  experiments,  it  was 
also  determined  that  the  trolley  angle  is  not  affected  by  the  structural  coupling  effects,  in 
contrast  with  the  trolley  twist  which  is  substantially  affected  by  structural  coupling. 


B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  Side  Trolley  Placement 

We  have  conducted  the  analysis  for  stern  placement  of  the  new  trolley-rail 
system.  However  there  are  cases  where  the  trolley-rail  system  may  be  placed  on  the  side 
of  the  ship  for  unloading  operations.  In  this  scenario,  the  hydrodynamic  coupling  effects 
of  the  ship  and  the  RRDF  will  be  significantly  different.  A  similar  analysis  should  be 
conducted  for  the  side  trolley  placement  to  detennine  the  relative  impact  of  the  structural 
coupling  effects. 

2.  Structural  Coupling  for  Trolley  Twist 

Given  the  knowledge  that  the  trolley  twist  is  substantially  affected  by  structural 
coupling  between  the  trolley  and  RRDF,  we  can  now  proceed  to  analyze  the  structural 
coupling  effects  utilizing  the  mathematical  model  developed,  so  that  appropriate 
measures  can  be  taken  to  ensure  that  cargo  transfer  can  be  conducted  at  sea  up  to  Sea 
state  5. 
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APPENDIX  A 


MATLAB  CODE 

%  Motions  calculation  for  CapeD/RRDF/trolley  interface. 

%  Long-crested  Bretchneider  seas. 

% 

i_seaway=input('Enter  0  for  PM  or  1  for  Bretchneider  spectrum  =  '); 
if  i_seaway  ==  1 

T  in  =input(’Modal  Period  (sec)  =  '); 
omega_m=2*pi/T_m; 
end 

i_stiffness=input('Enter  3  for  heave  stiffness,  4  for  roll,  or  5  for  pitch  =  '); 
k_stiffness=input('Enter  appropriate  normalized  stiffness  =  '); 

% 

x_ann_C  APED_  1  =-3  3  0 ; 
x_ann_C  APED_2=-3  3  0 ; 
y_ann_C  APED_  1  =- 1 0 ; 
y_ann_CAPED_2=+ 1 0; 
x_ann_RRDF_  1=20; 
x_ann_RRDF_2=2  0 ; 
y_ann_RRDF_  1  =- 1 0 ; 
y_ann_RRDF_2=+ 1 0; 
trolley_length=  1 00; 

% 

%  Load  processed  WAMIT  data  file 

% 

load  CAPED 
NH=25; 
i_NH=15; 
beta_incr=i_NH; 
imag_unit=sqrt(- 1 ); 

% 
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%  Set  the  frequencies  vector  (0.3  to  2.5  rad/sec) 

% 

NF=(size(frequencies)); 

NF=NF(1,1); 

w=frequencies; 

% 

%  Set  the  headings  vector  (0  to  360  degrees  in  i_NH  deg.  increments) 

% 

for  i=l:NH, 

heading(i)=i_NFI*(i-l); 

end 

% 

%  Set  added  mass  and  damping  matrices  and  forcing  vector 

% 

for  i=l:NF, 

i_string=num2str(i); 

A(:,:,i)=eval(strcat('admassfreq’,i_string)); 
B(:,:,i)=eval(strcat('addumpingfreq’,i_string)); 
for  j=l:NH, 
j_string=num2str(j); 
if  j==NH 

j_string=num2str(  1 ); 
end 

F(:,i,j)=eval(strcat('forcfreq’,i_string,’head’,j_string)); 

end 

end 

% 

%  Random  wave  calculations 

%  PM  spectrum  -  weighted  by  the  spreading  function 

% 

%  Loop  on  significant  wave  height 
iHS=0; 

for  HS=0. 5:0. 5:30, 
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HS; 

iHS=iHS+l; 
if  i_seaway  ==0 

omega_m=0.4*sqrt(32.2/HS); 

end 

A_s=(  1 .25/4)*(omega_mA4)*(HSA2); 

B_s=  1 .25  *omega_mA4; 

S_main=(A_s./w.A5).*exp(-B_s./w.A4); 

% 

%  Loop  on  sea  direction 

% 

ibeta=0; 

for  beta=0:beta_incr:360, 
beta_set=beta; 
ibeta=ibeta+l; 

S=S_main; 

heading_single=heading(ibeta); 

% 

%  Frequency  domain  response  is  given  by  x=  [-wA2*(A+mass)+i*w*B+C]A(-l)*F 
%  Zeroth  order  calculation  (no  trolley/RRDF  structural  coupling) 

% 

for  i=l:NF, 
w_single=w(i); 

A_bar=-w_single*w_single*(A(:,:,i))+hnag_unit*w_single!|sB(:,:,i)+C; 

F_bar=F(:,i,ibeta); 

x=inv(A_bar)*F_bar; 

% 

%  Ship  motions 

% 

surgeO_CAPED(i)  =  x(l); 
swayO_CAPED(i)  =  x(2); 
heaveO_CAPED(i)  =  x(3); 
rollOCAPED(i)  =  x(4); 
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pitchO_CAPED(i)  =  x(5); 
yawO_CAPED(i)  =  x(6); 

% 

%  RRDF  motions 

% 

surgeO_RRDF(i)  =  x(7); 
swayO_RRDF(i)  =  x(8); 
heaveORRDF(i)  =  x(9); 
rollORRDF(i)  =x(10); 
pitchORRDF(i)  =  x(l  1); 
yawO_RRDF  (i)  =  x(  1 2); 

% 

%  Motions/  at  trolley  end  -  GAPED  side 

% 

mO_t_C  APED_v_  1  (i)=heaveO_C  APED(i)- 
x_ann_C  APED_  1  *pitchO_C  APED(i). . . 

-yannCAPEDl  *rollO_CAPED(i); 

mO_t_CAPED_v_2(i)=heaveO_CAPED(i)- 
x_ann_C  APED2  *pitchO_C  APED(i). . . 

-y_ann_CAPED_2*rollO_CAPED(i); 

mO_t_CAPED_h_l(i)=abs(swayO_CAPED(i))+x_arm_CAPED_l*abs(yawO_CAPED(i)) 


mO_t_CAPED_h_2(i)=abs(swayO_CAPED(i))+x_arm_CAPED_2*abs(yawO_CAPED(i)) 


% 

%  Motions  at  trolley  end  -  RRDF  side 

% 

mO_t_RRDF_v_  1  (i)=heaveO_RRDF(i)-x_ann_RRDF_  1  *pitch()_RRDF(i). . . 
-y_arm_RRDF_  1  *rollO_RRDF  (i); 

mO_t_RRDF_v_2(i)=heaveO_RRDF(i)-x_ann_RRDF_2*pitchO_RRDF(i)... 
-y_ann_RRDF_2  *rollO_RRDF(i); 

rnO_t_RRDF_h_l(i)=abs(swayO_RRDF(i))+x_arm_RRDF_l*abs(yawO_RRDF(i)); 
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inO_t_RRDF_h_2(i)=abs(swayO_RRDF(i))+x_arm_RRDF_2*abs(yawO_RRDF(i)); 

% 

%  Trolley  relative  vertical  angular  displacement 

% 

mO_trolley_angle_v_l(i)=(abs(heaveO_CAPED(i))- 
xannC  APED_  1  *  abs(pitchO_C  APED(i)) . . . 

-y_ann_CAPED_l*abs(rollO_CAPED(i))-(abs(heaveO_RRDF(i))- 
x_ann_RRDF_  1  *  abs(pitchO_RRDF(i)) . . . 

-y_arrn_RRE)F_l*abs(rollO_RRDF(i))))/trolley_length; 

mO_trolley_angle_v_2(i)=(abs(heaveO_CAPED(i))- 
x_ann_C  APED2  *  abs(pitchO_C  APED(i)) . . . 

-y_ann_C  APED_2  *  abs(rollO_C  APED(i))-(abs(heaveO_RRDF(i))- 
x_ann_RRDF_2  *  abs(pitchO_RRDF(i)) . . . 

-y_arm_RRDF_2  *  abs(rollO_RRDF(i))))/trolley_length; 

% 

%  Trolley  relative  horizontal  displacement 

% 

mO_trolley_distance_h_  1  (i)=(mO_t_C  APED_h_  1  (i)-mO_t_RRDF_h_  1  (i)); 
mO_trolley_distance_h_2(i)=(mO_t_CAPED_h_2(i)-mO_t_RRDF_h_2(i)); 

% 

%  Calculate  transmission  force  due  to  trolley  end-points  motion  differential 

% 

if  i_stiffness==3 

f_diff(i)=C(9,9)*k_stiffness*(mO_t_RRDF_v_l(i)-mO_t_RRDF_v_2(i)); 

heave_trolley(i)=f_diff(i); 

roll_troliey(i)=0; 

pitch_troliey(i)=0; 

end 

if  i_stiffness==4 

f_diff(i)=C(10,10)*k_stiffness*abs((m0_trolley_angle_v_l(i)- 

m0_trolley_angle_v_2(i))); 

heave_trolley(i)=0 ; 

roll_trolley(i)=f_diff(i); 
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pitch_trolley(i)=0; 

end 

if  i_stiffness==5 

f_diff(i)=C(  1 1,1  l)*k_stiffness*abs((mO_trolley_angle_v_l(i)- 
m0_trolley_angle_v_2(i))); 

heave_trolley(i)=0 ; 

roll_trolley(i)=0; 

pitch_trolley(i)=f_diff(i); 

end 

end 

% 

%  First-order  correction 

%  Calculate  motions  again  incorporating  trolley  stiffness  correction 

% 

for  i=l:NF, 
w_single=w(i); 

A_bar=-w_single*w_single*(A(:,:,i))+imag_unit*w_single*B(:,:,i)+C; 

f_trolley=f_diff(i); 

m_trolley=f_diff(i)*20; 

F_bar=F(:,i,ibeta)+[0;0;0;0;0;0;0;0;heave_trolley(i);- 
rolltrolley(i);pitch_trolley(i );()]; 

x=inv(A_bar)*F_bar; 

% 

%  Ship  motions 

% 

surge_CAPED(i)  =  x(l); 
sway_CAPED(i)  =  x(2); 
heave_CAPED(i)  =  x(3); 
rollCAPED(i)  =  x(4); 
pitch_CAPED(i)  =  x(5); 
yaw_CAPED(i)  =  x(6); 

% 

%  RRDF  motions 
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% 

surge_RRDF(i)  =  x(7); 
sway_RRDF(i)  =  x(8); 
heaveRRDF(i)  =  x(9); 
rollRRDF(i)  =x(10); 
pitchRRDF(i)  =  x(l  1); 
yaw_RRDF(i)  =x(12); 

% 

%  Motions/  at  trolley  end  -  CAPED  side 

% 

m_t_CAPED_v_  1  (i)=heave_C  APED(i)-x_arm_CAPED_  1  *pitch_C  APED(i). . . 
-y_arm_C  APE  D_  1  *  ro  1  IGA  P  E  D(  i ) ; 

m_t_CAPED_v_2(i)=heave_CAPED(i)-x_arm_CAPED_2*pitch_CAPED(i)... 

-y_ann_CAPED_2*roll_CAPED(i); 

m_t_CAPED_h_l(i)=abs(sway_CAPED(i))+x_arm_CAPED_l*abs(yaw_CAPED(i)); 

m_t_CAPED_h_2(i)=abs(sway_CAPED(i))+x_arm_CAPED_2!|5abs(yaw_CAPED(i)); 

% 

%  Motions  at  trolley  end  -  RRDF  side 

% 

m_t_RRDF_v_l(i)=heave_RRDF(i)-x_ann_RRDF_l*pitch_RRDF(i)... 
-yarmRRDFl  *roll_RRDF(i); 

m_t_RRDF_v_2(i)=heave_RRDF(i)-x_ann_RRDF_2*pitch_RRDF(i)... 

-y_arm_RRDF_2*roll_RRDF(i); 

m_t_RRDF_h_l(i)=abs(sway_RRDF(i))+x_ann_RRDF_l*abs(yaw_RRDF(i)); 
m_t_RRDF_h_2(i)=abs(sway_RRDF(i))+x_ami_RRDF_2  *  abs(yaw_RRDF(i)); 
% 

%  Trolley  relative  vertical  angular  displacement 

% 

m_trolley_angle_v_  1  (i)=(abs(heave_CAPED(i))- 
x_ann_C  APED_  1  *  abs(pitch_C  APED(i)) . . . 

-y_ann_CAPED_l*abs(roll_CAPED(i))-(abs(heave_RRDF(i))- 
x_ann_RRDF_  1  *  abs(pitch_RRDF(i)) . . . 
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-y_arm_RRDF_l*abs(roll_RRDF(i))))/trolley_length; 

m_trolley_angle_v_2(i)=(abs(heave_CAPED(i))- 
x_ann_C  APED_2  *  abs(pitch_C  APED(i)) . . . 

-yannC  APED2  *  abs(roll_C  APED(i))-(abs(heave_RRDF(i))- 
x_ann_RRDF_2  *  abs(pitch_RRDF  (i)) . . . 

-y_arm_RRDF_2*abs(roll_RRDF(i))))/trolley_length; 

m_trolley_angle_v_twist(i)=m_trolley_angle_v_l(i)-m_trolley_angle_v_2(i); 

m_trolley_angle_v_average(i)=0 . 5  *  (m_trolley_angle_v_  1  (i)+m_trolley_angle_v_2(i)); 

% 

%  Trolley  relative  horizontal  displacement 

% 

m_trolley_distance_h_  1  (i)=(m_t_CAPED_h_  1  (i)-m_t_RRDF_h_  1  (i)); 
m_trolley_distance_h_2(i)=(m_t_CAPED_h_2(i)-m_t_RRDF_h_2(i)); 
end 

% 

%  Define  response  spectra 

% 

S_surge_CAPED  =  ((abs(surge_CAPED)).A2).*S'; 

S_surge_RRDF  =  ((abs(surge_RRDF)).A2).*S'; 

S_heave_CAPED  =  ((abs(heave_CAPED)).A2).*S'; 

S_heave_RRDF  =  ((abs(heave_RRDF)).A2).*S’; 

S_sway_CAPED  =  ((abs(sway_CAPED)).A2).*S’; 

S_sway_RRDF  =  ((abs(sway_RRDF)).A2).*S’; 

S_roll_CAPED  =  ((abs(roll_CAPED)).A2).*S’; 

S  roll  RRDF  =  ((abs(roll  RRDF)).A2).*S’; 

S_pitch_CAPED  =  ((abs(pitch_CAPED)).A2)AS’; 

S_pitch_RRDF  =  ((abs(pitch_RRDF)).A2).*S’; 

S_yaw_CAPED  =  ((abs(yaw_CAPED)).A2).*S’; 

SyawRRDF  =  ((abs(yaw_RRDF)).A2)AS'; 

S_m_t_CAPED_v_l  =  ((abs(m_t_CAPED_v_l)).A2).*S’; 

S_m_t_CAPED_v_2  =  ((abs(m_t_CAPED_v_2)).A2).*S’; 

S_m_t_C  APED_h_  1  =  ((abs(m_t_CAPED_h_l)).A2).*S’; 
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S_m_t_CAPED_h_2  =  ((abs(m_t_CAPED_h_2)).A2).*S’; 

S_m_t_RRDF_v_  1  =  ((abs(m_t_RRDF_v_l)).A2).*S'; 

S_m_t_RRDF_v_2  =  ((abs(m_t_RRDF_v_2)).A2).*S'; 

S_m_t_RRDF_h_  1  =  ((abs(m_t_RRDF_h_l)).A2).*S'; 

S_m_t_RRDF_h_2  =  ((abs(m_t_RRDF_h_2)).A2).*S'; 
S_m_trolley_angle_v_l  =  ((abs(m_trolley_angle_v_l)).A2).*S’; 
S_m_trolley_angle_v_2  =  ((abs(  mlro  1 1  cy_a  n  g  1  e_v_2 ) ) .  A2 ) .  *  S ' ; 
S_m_trolley_angle_v_twist  =  ((abs(m_trolley_angle_v_twist)).A2).*S'; 
S_m_trolley_angle_v_average  =  ((abs(m_trolley_angle_v_average)).A2).*S'; 
S_m_trolley_distance_h_  1  =  ((abs(m_trolley_distance_h_  1 )).  A2).  *  S'; 
S_m_trolley_distance_h_2  =  ((abs(m_trolley_distance_h_2)).A2).*S'; 

% 

%  Integral  initializations 

% 

SsurgeCAPEDi  =  0; 

SsurgeRRDFi  =  0; 

SheaveCAPEDi  =  0; 

SheaveRRDFi  =  0; 

S_sway_CAPED_i  =  0; 

SswayRRDFi  =  0; 

SrollCAPEDi  =0; 

SrollRRDFi  =0; 

S_pitch_CAPED_i  =  0; 

S_pitch_RRDF_i  =  0; 

S_yaw_CAPED_i  =0; 

S_yaw_RRDF_i  =  0; 

SnitC  APED_v_  1  _i  =  0; 

S_m_t_CAPED_h_l_i  =  0; 

S  m  t  RRDF  v_l_i  =0; 

SmtRRDFhli  =0; 

S_m_trolley_angle_v_l_i  =  0; 

Smtrolleydistancehli  =  0; 

S_m_t_CAPED_v_2_i  =  0; 
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S_m_t_CAPED_h_2_i  =  0; 

S_m_t_RRDF_v_2_i  =  0; 

S_m_t  RRDF_h_2  i  =0; 

S_m_trolley_angle_v_2_i  =  0; 

S_m_trolley_distance_h_2_i  =  0; 

S_m_trolley_angle_v_twist_i  =  0; 

S_m_trolley_angle_v_average_i  =  0; 

% 

%  Integral  S(w)*|RAO(w)|A2 

% 

for  i=2:l:NF, 

avg_value=0.5*(S_surge_CAPED(i)+S_surge_CAPED(i-l)); 

S_surge_CAPED_i  =  S_surge_CAPED_i+avg_value*abs(w(i-l)-w(i)); 

avg_valuc=0 .  5  *  (S_surge_RRDF  (i)+S_surge_RRDF(i- 1 )); 

S_surge_RRDF_i  =  S_surge_RRDF_i+avg_value*abs(w(i-l)-w(i)); 

avg_valuc=0 . 5  *  ( S_heave_C  APED(i)+S_heave_C  APED(i- 1 )) ; 

S_heave_CAPED_i  =  S_heave_CAPED_i+avg_value*abs(w(i-l)-w(i)); 

avg_value=0 . 5  *  (SheaveRRDF  (i)+S_heave_RRDF(i- 1 )); 

S_heave_RRDF_i  =  S_heave_RRF)F_i+avg_value*abs(w(i-l)-w(i)); 

avg_value=0 . 5  *  (S_sway_C  APED(i)+S_sway_C  APED(i- 1 )); 

S_sway_CAPED_i  =  S_sway_CAPED_i+avg_value*abs(w(i-l)-w(i)); 

avg_value=0 . 5  *  (S_sway_RRDF  (i)+S_sway_RRDF  (i- 1 )); 

S_sway_RRDF_i  =  S_sway_RRDF_i+avg_value*abs(w(i-l)-w(i)); 

avg_value=0.5*(S_roll_CAPED(i)+S_roll_CAPED(i-l)); 

SrollCAPEDi  =  S_rollCAPED_i+avg_valuc5|5abs(w(i- 1  )-w(i)); 

avg_value=0 . 5  *  (S_roll_RRDF(i)+S_roll_RRDF(i- 1 )); 

S_roll_RRDF_i  =  S_roll_RRDF_i+avg_valuc5|5abs(w(i-l  )-w(i)); 

avg_value=0.5*(S_pitch_CAPED(i)+S_pitch_CAPED(i-l)); 

S_pitch_CAPED_i  =  S_pitch_CAPED_i+avg_valueHtabs(w(i-l)-w(i)); 

avg_value=0 . 5  *  (S_pitch_RRDF(i)+S_pitch_RRDF(i- 1 )); 

S_pitch_RRDF_i  =  S_pitch_RRDF_i+avg_value*abs(w(i-l)-w(i)); 

avg_value=0.5*(S_yaw_CAPED(i)+S_yaw_CAPED(i-l)); 

S_yaw_CAPED_i  =  S_yaw_CAPED_i+avg_value*abs(w(i-l)-w(i)); 
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avg_value=0 . 5  *(S_yaw_RRDF(i)+S_yaw_RRDF(i- 1 )); 

S_yaw_RRDF_i  =  S_yaw_RRDF_i+avg_value*abs(w(i-l)-w(i)); 

% 

%  Calculations  for  point  "  1 " 

% 

avgval  uc=0 .  5  *  ( SnitCAP  ED_v_  1  (i)+S_m_t_CAPED_v_  1  (i- 1 )); 

SintC  APED_v_  l_i  =  SintC  APED_v_  1  _i+avg_value  *  abs(w(i- 1 )- w(i)); 
a  vg_  va  1  uc=() .  5  *  ( S_m_t_C  A  P  E  D_h_  1  (i)+S_m_t_CAPED_h_  1  (i- 1 )); 

SintC  APED_h_  l_i  =  SintC  APED_h_  1  _i+avg_value  *  abs(w(i- 1 )- w(i)); 
avg_value=0 . 5  *  (S_m_t_RRDF_v_  1  (i)+S_m_t_RRE)F_v_  1  (i- 1 )); 
S_m_t_RRDF_v_l_i  =  S_m_t_RRDF_v_l_i+avg_value!|5abs(w(i-l)-w(i)); 
avg_valuc=0 . 5  *  (S_m_t_RRDF_h_  1  (i)+S_m_t_RRDF_h_  1  (i- 1 )); 
S_m_t_RRDF_h_l_i  =  S_m_t_RRDF_h_l_i+avg_value5|5abs(w(i-l)-w(i)); 
avg_value=0 . 5  *  (S_m_trolley_angle_v_  1  (i)+S_m_trolley_angle_v_  1  (i- 1 )); 
S_m_trolley_angle_v_l_i  =  S_m_trolley_angle_v_l_i+... 
avg_valueH5abs(w(i- l)-w(i)); 

a  vg_  va  1  ue=0 . 5 5|!  ( Sintro  1 1  cycl  i  stan  c  c_h_  1  (i)+S_m_trolley_distance_h_  1  (i- 1 )); 
Smtrolleydistancehli  =  S_m_trolley_distance_h_l_i+... 
avg_valueH5abs(w(i-  l)-w(i)); 

% 

%  Calculations  for  point  "2" 

% 

avg_value=0 . 5  *(S_m_t_CAPED_v_2(i)+S_m_t_CAPED_v_2(i- 1 )); 
S_m_t_CAPED_v_2_i  =  S_m_t_CAPED_v_2_i+avg_value*abs(w(i-l)-w(i)); 
avg_valuc=0.55|!(S_in_t_CAPED_h_2(i)+S_m_t_CAPED_h_2(i- 1 )); 
S_m_t_CAPED_h_2_i  =  S_m_t_CAPED_h_2_i+avg_value*abs(w(i-l)-w(i)); 
avg_value=0 . 5  *  (S_m_t_RRDF_v_2(i)+S_m_t_RRDF_v_2(i- 1 )); 
S_m_t_RRDF_v_2_i  =  S_m_t_RRDF_v_2_i+avg_value*abs(w(i-l)-w(i)); 
avg_valuc=0 . 5  *(S_m_t_RRDF_h_2(i)+S_m_t_RRDF_h_2(i- 1 )); 
S_m_t_RRDF_h_2_i  =  S_m_t_RRDF_h_2_i+avg_value*abs(w(i-l)-w(i)); 
avg_value=0 . 5  *  (S_m_trolley_angle_v_2(i)+S_m_trolley_angle_v_2(i- 1 )); 
S_m_trolley_angle_v_2_i  =  S_m_trolley_angle_v_2_i+... 
avg_value*abs(w(i-l)-w(i)); 
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avg_value=0.5*(S_m_trolley_distance_h_2(i)+S_m_trolley_distance_h_2(i-l)); 
S_m_trolley_distance_h_2_i  =  S_m_trolley_distance_h_2_i+... 
avg_valuc*abs(w(i- 1  )-w(i)); 

% 

avg_value=0 . 5  *  (S_m_trolley_angle_v_twist(i)+S_m_trolley_angle_v_twist(i- 1 )); 
S_m_trolley_angle_v_twist_i  =  S_m_trolley_angle_v_twist_i+... 
avg_value*abs(w(i-l)-w(i)); 

avg_value=0.5*(S_m_trolley_angle_v_average(i)+S_m_trolley_angle_v_average(i-l)); 
S_m_trolley_angle_v_average_i  =  S_m_trolley_angle_v_average_i+. . . 
avg_value*abs(w(i-  l)-w(i)); 
end 
% 

%  RMS  values 

% 

RMS_surge_CAPED(ibeta,iHS)=  sqrt(S_surge_CAPED_i); 
RMS_surge_RRDF(ibeta,iHS)  =  sqrt(S_surge_RRDF_i); 
RMS_heave_CAPED(ibeta,iHS)=  sqrt(S_heave_CAPED_i); 
RMS_heave_RRDF(ibeta,iHS)  =  sqrt(S_heave_RRDF_i); 
RMS_sway_CAPED(ibeta,iHS)  =  sqrt(S_sway_CAPED_i); 
RMS_sway_RRDF(ibeta,iHS)  =  sqrt(S_sway_RRDF_i); 
RMS_roll_CAPED(ibeta,iHS)  =  sqrt(S_roll_CAPED_i); 
RMS_roll_RRDF(ibeta,iHS)  =  sqrt(S_roll_RRDF_i); 
RMS_pitch_CAPED(ibeta,iHS)=  sqrt(S_pitch_CAPED_i); 
RMS_pitch_RRDF(ibeta,iHS)  =  sqrt(S_pitch_RRDF_i); 
RMS_yaw_CAPED(ibeta,iHS)  =  sqrt(S_yaw_CAPED_i); 
RMS_yaw_RRDF(ibeta,iHS)  =  sqrt(S_yaw_RRDF_i); 
RMS_m_t_CAPED_v_l(ibeta,iHS)  =  sqrt(S_m_t_CAPED_v_l_i); 
RMS_m_t_CAPED_h_l  (ibeta,iHS)  =  sqrt(S_m_t_CAPED_h_l_i); 
RMS_m_t_RRDF_v_l (ibeta,iHS)  =  sqrt(S  m_t  RRDF  v  l  i); 
RMS_m_t_RRDF_h_l(ibeta,iHS)  =  sqrt(S_m_t_RRDF_h_l_i); 
RMS_m_trolley_angle_v_l(ibeta,iHS)  =  sqrt(S_m_trolley_angle_v_l_i); 
RMS_m_trolley_distance_h_l(ibeta,iHS)  =  sqrt(S_m_trolley_distance_h_l_i); 
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RMS_m_t_CAPED_v_2(ibeta,iHS)  =  sqrt(S_m_t_CAPED_v_2_i); 

RMS_m_t_CAPED_h_2(ibeta,iHS)  =  sqrt(S_m_t_CAPED_h_2_i); 

RMS_m_t_RRDF_v_2(ibeta,iHS)  =  sqrt(S_m_t_RRDF_v_2_i); 

RMS_m_t_RRDF_h_2(ibeta,iHS)  =  sqrt(S_m_t_RRDF_h_2_i); 

RMS_m_trolley_angle_v_2(ibeta,iHS)  =  sqrt(S_m_trolley_angle_v_2_i); 

RMS_m_trolley_distance_h_2(ibeta,iHS)  =  sqrt(S_m_trolley_distance_h_2_i); 

RMS_m_trolley_angle_v_twist(ibeta,iHS)  =  sqrt(S_m_trolley_angle_v_twist_i); 

RMS_m_trolley_angle_v_average(ibeta,iHS) 

sqrt(S_m_trolley_angle_v_average_i); 

end 

end 

% 

figure(l) 

[th,r]=meshgrid((0:beta_incr:360)*pi/ 1 80,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.0:0. 1:2.0]; 

contour(X',Y',RMS_m_trolley_angle_v_average*  1 80/pi, c_p),caxis([0  2]),colorbar 
title(’RMS  value  of  Trolley  Angle  (deg)') 

% 

figure(2) 

[th,r]=meshgrid((0:beta_incr:360)*pi/l  80, 0.5:0. 5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.0:0. 1:2.0]; 

contour(X',Y',RMS_m_trolley_angle_v_twist*  1 80/pi, c_p),caxis([0  2]),colorbar 
title(’RMS  value  of  Trolley  Twist  (deg)') 
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